Reply to "Comment on 'Stimulated Raman adiabatic passage from an atomic to a 

molecular Bose-Einstein condensate'" 

P. D. Drummond and K. V. Kheruntsyan 
ARC Centre of Excellence for Quantum-Atom Optics, Department of Physics, 
University of Queensland, Brisbane, Qld 4072, Australia 

D. J. Heinzen and R. H. Wynar 
Department of Physics, University of Texas, Austin, Texas 78712 
(Dated: February 2, 2008) 

In the Comment by M. Mackie et al. [arXiv: physics/0212111 v. 4], the authors suggest that the 
molecular conversion efficiency in atom-molecule STIRAP can be improved by lowering the initial 
atomic density, which in turn requires longer pulse durations to maintain adiabaticity. Apart from 
the fact that the mean-field approximation becomes questionable at low densities, we point out that 
a low-density strategy with longer pulses has several problems. It generally requires higher pulse 
energies, and increases radiative losses. We also show that even within the approximations used 
in the Comment, their example leads to no efficiency improvement compared to our high-density 
case. In a more careful analysis including radiative losses neglected in the Comment, the proposed 
strategy gives almost no conversion owing to the longer pulse durations required. 
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The authors of Ref. jl| revisit our earlier work on 
stimulated Raman adiabatic passage (STIRAP) from 
an atomic to a molecular Bose-Einstein condensate [2|. 
They suggest that the disruptive role of dephasing due 
to particle-particle interactions can be reduced by lower- 
ing the initial atomic condensate density. This strategy 
requires a simultaneous increase of the duration of the 
Raman pulses. As a result the authors claim to obtain 
improved conversion efficiency, which reaches about 51% 
for the optimized density. This is about 5% larger than 
the highest efficiency of ^46% found in Ref. [2( , and is 
achieved by using a 100 times lower atomic density and 
about 100 times longer pulse duration. 

Apart from the fact that mean-field theory becomes 
questionable for describing low density condensates, we 
start by emphasizing that any comparison of efficiency is 
only meaningful within a well-defined situation of equal 
atom numbers and laser powers (or pulse energies) Q. 
This is because the optimum efficiency is a strong func- 
tion of maximum Rabi frequency. Given these con- 
straints, together with increased radiative losses at longer 
pulse durations, we show here that there are three signif- 
icant problems with the results of Ref. 0] . 

Laser power. The first problem is that a low density 
strategy requires a higher laser power to obtain an equal 
laser intensity, since the laser waist size must be enlarged. 
This is due not only to the larger size of the condensate, 
for fixed atom number, but also to the fact that the ex- 
periment becomes more sensitive to the inhomogeneity 
in the AC Stark shifts from the photoassociation laser 
beams, which set a practical lower limit on the beam 
waists. If such high power lasers to implement the low 
density strategy of Ref. flj are available, they can sim- 
ply be utilized to increase the maximum Rabi frequency 
in the original proposal of Ref. @ - which greatly in- 
creases the conversion efficiency. Using fewer atoms, as 



suggested in the Comment, invalidates the comparison. 

Parameter normalization and optimization. The next 
problem is that equations (la)-(lc) of Ref. pj and the 
parameters used do not correspond to those of Ref. • 
The difference is in the starting value of the parameter 
Xo an d the fact that the respective equations of Ref. @ 
contain an extra factor 1/V2 in front of this coupling 
term. Hence, the comparison is made between different 
Rabi frequencies, not just a lower atomic density. 

Thus, the results of the Comment should be compared 
with uniform trap results of Ref. Q obtained with \/2 
times larger Rabi frequency fli (see Eqs. (6) of Ref. 
0). In this case, simulation of the respective equations 
of Ref. with a factor of y/2 larger peak value of fii, 
zero two-photon detuning, and the same original values 
for the remaining parameters (i.e. higher atomic density 
po = 4.3x 10 14 cm -3 , apulse duration of T — 0.2 ms, and 
a 50 times larger peak value of O2 than fii) gives 67% 
conversion efficiency. This is higher than the maximum 
efficiency of 51% found in the Comment pj. 

In other words, the claimed improvement in conversion 
is obtained not just by altering the initial atomic density 
and pulse duration, but by changing other parameter val- 
ues as well. When directly comparable Rabi frequencies 
are used, the relative efficiency is in fact worse in the 
suggested low-density strategy. 

We also point out that the conversion efficiency of 16% 
seen in Fig. 1 of the Comment at po = 4.3 x 10 14 cm -3 
should not be confused with our result of 67% at the 
same density. The reason for this seeming disagreement 
is that the authors of the Comment use essentially a one- 
parameter (density) optimization. This means that other 
parameter values like pulse duration are not necessarily 
optimized at all densities. For example, the pulse du- 
ration is set to T = 5 x 10 3 /xo in a ll cases, which for 
po = 4.3 x 10 14 cm -3 gives T = 2.4 ms. In contrast to 
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this, our original optimization procedure is carried out 
with respect to three parameters - the pulse duration, 
pulse offset, and the two-photon detuning. The 67% effi- 
ciency obtained in our case uses an optimum pulse dura- 
tion of T = 0.2 ms, which is much shorter than used in 
the Comment at this density Q. 

It should be noted that 67% efficiency is higher than 
the maximum of 46% we found previously |2|. This is 
simply due to the improved adiabaticity following ob- 
tained at higher Rabi frequency. 

Radiative losses. We now address the most detrimen- 
tal outcome of the longer pulse durations employed in the 
Comment. As the timescale for STIRAP is increased, in- 
coherent radiative losses become increasingly important 
due to additional radiative couplings that are ignored in 
the model of Ref. 0. We show here that once included 
in the model, they can lead to dramatic losses of the 
ground-state molecules. Our results for typical param- 
eter values show almost zero final population. Thus, in 
the suggested case of 87 Rb, the proposed strategy results 
in a dramatic reduction rather than improvement in the 
conversion efficiency. 

There are two types of additional radiative couplings 
that are neglected in the Comment: (1) couplings that 
occur even within the simple three-level model, and (2) 
couplings to other vibrational levels in the excited molec- 
ular potential. 

The couplings within the three-level model included 
in our original paper 0, @ are due to non-resonant in- 
teractions of the two Raman lasers. These lead to the 
following, more complete equations: 
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the atoms, excited molecules, and stable molecules in the 
ground potential, respectively. 

For the benefits of the readers, we use the same no- 
tations as in Ref. Here, the density-dependent cou- 
pling x describes free-bound transitions between atom 
pairs and excited molecules due to the Raman laser at 
frequency uj\, which is the second pulse in the coun- 
terintuitive STIRAP sequence. The Rabi frequency fl 
describes the bound-bound transitions between excited 
and ground-state molecules due the first Raman pulse 
at frequency uj 2 . In addition to these couplings, which 
are the primary transitions in any STIRAP scheme, the 
above equations include the coupling of atoms to excited 
molecules due to the ix>2-laser, and similarly the coupling 
of excited and ground-state molecules due to the wi-laser. 
The respective coupling constants are x an d fl, and in 
the case of A = 5 = these non-primary or 'partner' 
transitions are detuned by U12. 

Explicit parameter values here are taken as in Ref. 
to make the comparison valid: A = 0, 5 = xo, 7s = 
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s 1 , and the initial atomic density is p = 4.3 x 10 12 cm 3 
so that xo 



2.1 xlO 5 



T = 5xl0 3 /Xo 



s '.In addition, fl = 50xo = 10' 
24 ms, Di = 4.5T, and D 2 = 2.5T. 

and 



Finally, we take xo = 10 4 s - \ fl Q = 2.3 x 10 s s~ 
u>i2 = —5.1 x 10 10 s _1 which are close to the typical 
calculated values corresponding to the spectroscopically 
most favorable case treated in Ref. Q- The values of 
the new parameters xo and CIq are chosen in a favorable 
manner [6j , and our conclusions would remain valid if the 
above mentioned 'missing' factors of l/v2 were restored 
self-consistently. 

Using the above parameters and simulating Eqs. 
, gives a final population of the ground-state molecules 
of \g\ 2 ~ 0, while the peak value during the pulse se- 
quence reaches only \g\ 2 ~ 2.5 x 10 ~ 6 . This implies es- 
sentially zero conversion efficiency. 

The reason for this dramatic result is that the newly 
formed ground-state molecules arc being still illuminated 
by the laser u>i during the second Raman pulse. As a 
result, they experience radiative losses at a rate of 
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X (t) =Xoexp [-(t-Drf/T 2 ] , 

X(t) =X cxp [-{t-D 2 ) 2 /T 2 } , 

n(i) = ft exp [~(t - D 2 ) 2 /T 2 } , 

fi(t) = O exp [-(i - D x ) 2 /T 2 } , 

while CJ12 = u>i — LU2 is the frequency difference between 
the two Raman lasers, A (5) is the two-photon (inter- 
mediate) detuning, and Aj„- represent atom-atom, atom- 
molecule, and molecule-molecule s-wave scattering inter- 
actions. The complex amplitudes a, 6, and g represent 
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Thus, the characteristic time scale for losses is 1/T e ff ~ 
2.6 ms, which is much smaller than the pulse duration 
T ~ 24 ms employed in Ref. [lj. These radiative losses 
are negligible for the much shorter (sub-millisecond) 
pulses treated in Ref. In this case, r e //T <C 1 and 
the role of the non-primary transitions is negligible. 

One might argue that targeting lower-lying vibrational 
levels in the ground molecular potential would give larger 
detuning W12, thus making T e ff smaller. This approach, 
however, suffers from the fact the respective bound- 
bound Franck-Condon overlap integrals typically become 
smaller and hence give even smaller values of the Rabi 
frequency Qq. In addition, a detailed multilevel analysis 
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reveals that these lower lying levels do not correspond to 
the most favorable case, once we take into account the 
entire set of necessary conditions for efficient conversion 
(see Eqs. (48)- (54) and the typical parameter values in 
Table IV of Ref. [2j ) . The reason for this is that increas- 
ing the detuning u>x2 will eventually bring the frequency 
of one of the lasers (or both) to a nearby resonance in the 
excited potential, thus giving rise to additional induced 
losses just as in Eq. Q}. 

The induced molecular losses due to the couplings to 
other vibrational levels in the excited potential can be 
modeled by the following additional term in the right- 
hand side of Eq. 

ig =(■■■)- 2 F 25, (5) 

where 

r 2 = r « e - 2 (*-^) 2 /T 2 + r (2) e - 2(t -D 3 )Vr^ (6) 

and the expressions for are given by Eq. (43) of Ref. 

. Here, the most disruptive loss coefficient is which 
is due to the bound-bound transitions during the sec- 
ond Raman pulse (laser wi). We take — 400 s _1 
here, which is more favorable than the value obtained in 
Ref. 0. This gives a characteristic time scale for losses 

l/F^ ~ 2.5 ms, which is much shorter than the pulse 
duration T ~ 24 ms. Not surprisingly, the simulation of 
Eqs. CEJ-© with the additional loss term, Eq. (JSJ, and 
with xo an d f^o set to zero, gives again almost zero final 
population of the ground-state molecules. 

An examination of the conditions stated in Eqs. (48)- 
(54) of Ref. and of the characteristic values of the 
respective coefficients in Table IV reveals that these re- 
sults could have been expected. The entries in lines 5 and 
12 of Table IV give the typical values of and T e ff. 
With pulse durations of T ~ 24 ms, these coefficients do 
not satisfy the conditions T -C I and T e ffT -C 1, and 



therefore the induced losses are not negligible. Similarly, 
even the induced atomic loss coefficient = 51 s" 1 
gives a^T ~ 1.2 (with T — 24 s _1 ) and hence it cannot 
be neglected either. 

We should note that all of the calculations we pre- 
sented have been carried out within a simplified model 
of the molecular structure, which assumes an 87 Rb2 
molecule with a single ground and single excited elec- 
tronic state. We include vibrational, but neglect rota- 
tional and hyperfine structure. As such, our conclusions 
- while correct qualitatively - can be further improved 
for quantitative purposes. Nevertheless, in contrast to 
the calculations of Ref. 1], ours do include the impor- 
tant physical effects in a reasonably realistic way, within 
the limitations of mean-field theory. We feel that the gen- 
eral conclusions we have presented, such as the poor scal- 
ing of the molecular conversion efficiency with decreasing 
density, and the lack of very substantial improvement in 
conversion with increasing molecular binding energy, arc 
likely to be present in any STIRAP experiment with al- 
kali dimers. 

To summarize, the photoassociation strategy employed 
in Ref. .1] requires higher laser powers and pulse energies. 
Even neglecting losses, it does not lead to improvements 
in efficiency relative to the original proposal, provided 
the maximum Rabi frequencies are kept equal. It also 
fails to capture the physics of induced radiative losses 
relevant for long pulse durations. The proposal of using 
lower atomic densities and longer pulses for achieving 
higher conversion efficiencies in atom-molecule STIRAP 
gives in fact the opposite effect, once these additional 
loss channels are taken into account. We find that with 
typical parameter values, a more complete model results 
in almost zero conversion, under the proposed conditions. 
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[3] The Rabi frequencies chosen in the original proposal 
would correspond approximately to a 100 (im waist and 
1 W power per Raman channel, typical of current cw 
Ti:sapphire laser beams. If the beam waist size is much 
smaller than this, then its AC Stark shift becomes very 
inhomogeneous, which creates new problems. 

[4] It is possible that a better optimization procedure might 
work better at lower densities, within the model employed 
in the Comment. However, we have not investigated this in 
detail, as the low density strategy will have to necessarily 
employ longer pulse durations to maintain adiabaticity, in 
which case the model itself becomes physically unrealistic 



due to the neglected induced radiative losses. 

[5] Eq. (29) of Ref. contains a typographical error. Here, 
the exponential term exp(— iwi2i) in the first line should 
read as exp(iuji2t) . The same term in Eq. (57) is correct. 
In addition, in the caption to Table III, |ii,3| =0.1 should 
be replaced by |Xa,3 1 = 0.1, and the quoted values of 8 
should refer to —8 instead. The value of |/2,3| is correct 
elsewhere in the text, and these typographical errors do 
not affect the results presented in Ref. 0. 

[6] On evaluating the values of the parameters xo and f2o , one 
has to pay attention to the fact that these have to be cal- 
culated in conjunction with xo and Qq and that they can 
not be varied independently. The reason is that these cou- 
plings originate from the values of the bare electronic Rabi 
frequencies for the two transitions, f2$ eI '°' ) and O.^ l ' \ and 
the free-bound and bound-bound Franck-Condon over- 



lap integrals, ii,3 and 72,3, using the notations of Ref. xo <x ^2 ' -^1,3 an d cx 0.[ el '°^ h,3- 

0- Accordingly, xo oc Q,^ 1 ' ^ Ii^, Q.o oc f2^ e!,0 'i"2,3, while 



